Fluctuations of Random Noise Power* 

By D. SLEPIAN 

{Manuscript received September 4, 1957) 

The probability distribution of the power, y, of a sample of Gaussian 
noise of time duration T is considered. Some general theory is presented 
along with curves for the cumulative distribution and probability density of 
y for several different power spectra and values of T. 

I. INTRODUCTION 

A random quantity of interest in many communication and detection 
systems is the average power, 

fin 

y = U ' N\t) dt, (1) 

J J-T/2 

of a sample of finite time duration, T, of a Gaussian noise, A''(/). This 
(|uantity has lieen discussed in some detail by Rice in his classic paper 
where he obtains expressions for the first few moments of y and an ap- 
proximate probability density function. 

In tliis paper the exact probal)ility density function, f{y), and the 
cumulative distribution function, F{y), of the average power are com- 
puted for a nimiber of ergodic Gaussian noises and for a number of 
values of T. The results are presented as a series of curves which are dis- 
cussed in the next section. It is hoped that they will be of use to tho.se 
designing specific systems. 

II. SUMMARY OF COMPUTATIONAL RESULTS 

V\g. 1 shows the pi'oi]ability density function, fiy), for the random 
varialile // of ec[Uatiou (1) when .V(0 has mean zero and power spectrum 

2a 
a- + 'iir-f- 

Noise with this spectrum will he referred to as RC noise (see 5.1). 

* The research rei)orte(l here was .supported in part bj- the Office of Naval 
Research under contract Nonr 210(00). 
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Fig. 1 — Probability density, /(</), foi- R(: noise. 

The curves are labelled by values of /3 = aT j'l. The curve marked 
/3 = is the probability density function for \j = N'(l). Fig. 2 shows the 
corresponding cumulative distribution functions, Fiy). 

For any > 0, as y approaches zero, /(//) and F{y) approach zero 
more rapidly than any power of y. 

As becomes large, the density function /(?/) peaks up around unity 
which is the average power of N{t). The variance of y is given by 
(2;8)~^[4^ — 1 + e~^\ It approaches zero for large /5 like ^ \ 

Figs. 3, 4 and 5 show f{y) when N{t) has mean zero and power spec- 
trum 






Wo 



+ 



— ) {w — Wd ; 



(2) 



O) = 2x/, - ^ ^ / ^ 



Noise with this spectrum will be referred to as RLC Noise (see 5.2). 
The figures are respectively for the cases Q = 1, 10 and 100. The curves 
are labelled by values of s - wqT. The curves marked s = are the 
density function for y = N^{t). The corresponding cumulative density 
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Fig. 2 — C'umuhilivc disli'ihution, /''(//), for KC noiae. 



functions, Z'Xy), i^i'e Khmvn on Figs. 6, 7 and 8. The spectra for Q = I, 
10 and 100 are plotted on Fig. 9. 

Vov any s > and for any finite Q > 0, as y approaches zero, both 
/(/y) and F(y) approach zero more rapidly than any power of y. 

Vov any fixed Q, as s becomes large, the density function f(y) peaks 
up around unity which is the average power of A'^(0- The variance of y 
is given liy 



^' %'[-! + ^--^ + ^, ^^^' \ V4Q^], 



(3) 



s 



For fixed Q, it approaches zero for large s like 2Q/s. 

If, however, s = wiJ' is held fixed and Q is permitted to increase, Figs. 
'A, -i and o show that f{y) becomes less concentrated; that is, with 
fixed integration time and fixed resonant frequency, fluctuations in 
power become more pronounced as the relative width of the spectral 
peak is decreased. Indeed, one has 



lim <t' ^ 1 + 

Q-»ao 



sm s 



so that 
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lira liin a = 1, 



8-KlO Q-*CO 



whereas, as already noted, 



lim lim. o- = 0. 



In the limit Q = «:> , the Gaussian noise can be taken to be the single 
frequency ensemble N{i) = a cos mi + b sin uid, where a and h are in- 
dependent normal variates with mean zero and variance unity. The 
density for y in this ease is 



/(?y) = sec (pe 



—y Bv.a^'fi 



Jaiiy tan <p sec (p) 



where sin <f = sin i/s and Jo is the usual Bessel function (see Appendix 
1). This density is plotted for several values of s in Fig. 10. It is to be 
noted that this limiting noise, although stationary, is not ergodic. It is 
this fact that causes the variance of y to be bounded away from zero as 
g _^ oo . Quite generally, if N{1) has a purely continuous spectrum, the 
variance of y will approach zero as the integration time becomes infinite. 
If the spectrum of A^{/} has line components, this vnW not be the case. 

It is not difficult to give a qualitative argument as to why power fluc- 
tuations in a fixed time interval increase as the power spectrum becomes 



2.0 



f(y) 




Fig. 3 — Probability density, f(y), for RLC noise, Q = 1.0, s = uaT. 



FLUCTUATIONS OF RANDOM NOISE POWER 



167 



f(y) 




Fig. 4 — Probaljility deiiwitj', J[y), for RLC noise, Q = 10, s = moT. 
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Fig. 5— Probability density, /(y), for RLC noise, Q = 100, s = woT. 
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1.2 




Fig. 6 — Cunmliitive distribution, Fiy), for RLC noise, Q = 1.0, s = w^T. 

more peaked. Noise with the power spectrum (2) can be thought of as 
the noise voltage produced across the resistor in a series RLC circuit 
when the applied \'oltage to the circuit is white Gaussian noise. The 
larger the Q of the circuit, the more it tends to "ring" in response to an 
impulse input; i.e., the longer the transients persist. An atypical excur- 
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Fig. 7 — Cumulative distribution, F(ij), for RLC noise, Q = 10, s = wuT. 
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Fig. 8 — Cumulative distriliutioii, /•'(]/), for RLC noii^c, Q — 100, s = uoT. 

sion of the input voltage will therefore have a longer lasting effect in the 
output of a circuit with a large Q than in a circuit with a small Q. To 
obtain the same variance in power, then, the integration time must be 
longer in the circuit with the large Q value. It would seem reasonable to 
expect this argument to apply for any peaked spectrinn, not solely for 
(2).* 

If the Q of the spectnmi (2) is increased, how much must the integra- 
tion time be increased to maintain roughly the same power fluctuations? 
From (3), it is seen that for large Q, a' is approximately 2t [r — 1 + e^^], 
i.e., a function of 

r = - = ^ T 

Q Q 

alone. Now Q measures the relative sharpness of the spectral peak, so 
that Wd/Q is a measure of the absolute width of the peak in radians/sec. 
As a rough rule, then, power measurements from different members of 
the family (2) will have the same fluctuations if their products "integra- 
tion time" times "absolute spectral bandwidth" are the same. Fig. 11 
shows ff^ as a function of r for Q = 1 , 10, and 100. That t is a good meas- 
ure of the fluctuation in power can also be seen by comparing the / 
curves of equal r value in Figs. 3, 4 and 5. They are almost identical. 



* It aoems to be very difficult to make any other qualitative statements re- 
garding the relation between the shape of the noise spectrum and the density 
function for y. 
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On Fig. 11 the variance of y for bandpass noise with spectrum 

I / ± /o I ^ 5 
to, I / ± /o I > fi 



w{j) = 



45' 



is plotted versus 



2TrU 



Here the Q of a bandpass circuit is defined by 



200 



O-oW 




3.5 



Fig. 9 — RLC spectra, Q = 1, 10, 100. 
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Fig. 10 — Probability den.sity, /(y), for RLC noise, (3 = =» . 



jind measures the relative width of the spectrum. This definition of Qi, 
ciiu.ses the a' curves of this noise power to agree asymptotically mth 
those of the RLC noise power; namely a^ '^ 2/r in both cases. Again, 
when it is not too small, r seems to be a good measure of power fluctua- 
tions. The variance in this bandpass case is given by 



.r = 4 f (I - y) 
Jq 



cos Qbrysm — - 

Trry 
2 



dy 



which can be readily evaluated in terms of Si and Ci functions. The curve 
tor Qh = 100 coincides .so closely with the curve for Qt = 10 it could not 
l)e shown on Fig. 11. 

The asymptotic agreement of the variance of noise power for band- 
pass and HLC noise permitted defining the Q of the bandpass circuit as 
Qi> = 'r(/ii/25). These same considerations suggest defining the band- 
width ir of the RLC spectrum by W ^ (i}o/2Q. For, in the bandpass 
('as(^, T = 2(25)7' which is 2T times the bandwidth of thespectnun. For 
the RLC noise, r = con'/'/C^ = 2T(iOo/2Q), whence the definition of IT 
follows. 
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Fig. 11 — a- for RLC noise and bamlpass noise, 0=1, 10, 100. 



The curves shown in Figs. 1-8 are believed to be accurate to two sig- 
nificant figures. For comparison, some points computed from Rice's 
approximate formula for f(y) (e<iuatiou 3.9-20 of [Kef. 1]) are shown on 
Fig. 8. Rice's formula is seen to fit the tails of /{//) well for large ij, but 
the central portion of the distril)ution is given accurately only for large 
values of t. However, the approximate cumulative distribution obtained 
by integrating Rice's formula agrees quite well with F{y) for a wide 
range of t values as is seen in Fig. G. 
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The approximiition in tiueHtion assumes a x^ type distribution 

(7,/2)-l -(!//2g2) 

Ky) = '^ — ' 



(2g^)'"T '^ 



(!)■ 



The parameters g and n are chosen to make the first two moments of 
this density agree with the true first two moments of y. That is, for the 
normalization 7?y/ = 1 adopted here, the equations g^re = I and2ng* = a^ 
serve to determine fj and n. These formulae give n = Ija^. Since for 
large t, q^ '^ 11 t for bandpass noise, n ^^ r = 2(25)7". That is, for large 
T, the bandpass noise acts like a x" variate with 2(25)7' degrees of free- 
dom in agreement with an argument easily derived from the sampling 
theorem. 

III. GENERAL THEORY 

Let N{V) be a Gaussian noise with mean zero and covariance 

pit, i') = E\N{t)I^{i')\ 

where as usual 7-7 denotes expectation. In studying properties of iV(/) in 
a finite time interval, say ( — 7'/2, 7'/2), it is convenient to make an ex- 
pansion in terms of an orthonornial set of functions, <pJS)^ n = 0, 1, 2, ... . 
We write 



where 



and 






T 

2 




/(, - ( N{t)<p,{t) (It, 




i ^ 0,1,2, 






'■,i- 0,1,2, . 



J-T/2 



As is well known,' it is particularly convenient in this description of 
the noise to choose as the orthonormal set, tpi , the solutions of the 
homogeneous Fredholni e(iuation with pit, i') as kernel. That is, the ^'s 
are chosen so tiiat 

x,^,(0 - f pC i')<pM') <n', I /. I ^%, i = 0,1,2, • ■ ■ . (4) 

J-Tl-l 2 
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For, with this choice of the ^'s, it is easily shown that the Ui are inde- 
pendent Gaussian variates with mean zero and variance Eiui ) = Xi , 
i = 0, \,2, ... . We assume in all that follows that the X's are so labelled 

that Xo ^ Xi ^ Xs ^ . . ■ ■ 

Consider now the average power, y, of a finite sample of the noise. It 

follows that 

y = )r N\t)dt = lj,T.n<' 

1 J-Tl2 i 



(5) 



= 2 "-^'^ 



where 



and 



Hi 



V^i 



X. 



(6) 



Equation (5) exhibits y as a linear combination of independent random 
variables. The Xi are independent Gaussian variables all with mean zero 
and variance unity. The characteristic function, C(^), for y then follows 
readily. One has 



(7) 

^ n (1 - 2^w-"^ 



j=0 



Here, as throughout this paper, the positive square root of a complex 
quantity is taken to have an angle between -(t/2) and +(t/2) radians 
(the cut line is along the negative real axis). 

From the characteristic function (7), the semi-invariants of y can be 
calculated. By definition^ of the semi-invariants, kj , 

1 C! 



From (7) and the expansion 



CO n 
X 



log (1 - a;) = -X) — . 

n=l ft' 
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it follows that 

log m = -5 i: log (1 - 2i£a,) ^Itt ^' 

V (^■^)' 

= Z- -77]- «^- 

where 

K, = (A:- 1)! 2*^-^1:7^ a/. (8) 

From the semi-invariants, the moments of y can be foimd as in Refer- 
ence 3. 

The formula (8) for the semi-invariants can be put in a convenient 
form not involving the aj exphcitly. From the well known expansion 

p(t, t') = E" X;^y(0^i(i') 
and the orthonormal properties of the ^'s, one finds 

i'' J- Tl2 

where the iterated kernel p (t, t') is defined by 

p''\t,i!) = p{t,t'), 
p'^'it, t') = f pit, x)p"'-'\x, t') dx, n - 2, 3, • • ■ . 

J-TiZ 

The determination of the higher order iterated kernels generally becomes 
difficult in practice. 

The expression (9) is of the form conjectured by Rice^ on the basis of 
computing the first four semi-invariants of y. The formula (7) was given 
by Kac and Siegert' and (9) was noted by Arthur'' in a special case in 
connection with the analysis of a frequency discriminator. 

The probability density function for y is obtained as the Fourier 
transform of C'{^), 

f(y) - ~ / ^ — (10) 

and the cumulative distribution function can be written as 

Fiy) = 1-1 fix) dx. in) 
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Much of the remainder of this paper will be concerned with evaluating 
(10) and (11) for specific noises. 

IV. COMPUTATIONAL FORM FOR fiy) 

The evaluation of the integral (10) presents many difficulties even 
with modern computing machinery. From the physical origins of the 
problem under discussion, it is clear that for smaU values of T, f(y) 
must be a rather broad function (non-localized), whereas for large values 
of T it must approach a 5-function centered at the point y = p(0, 0) 
when the noise is assumed ergodic. The behavior of (10) therefore de- 
pends in detail on the manner in which the a, approach zero with in- 
creasing j. 

One seemingly attractive approach to the problem is to truncate the 
sum in (5) at i = M and correspondingly obtain a product with ./ run- 
ning from to ]\f in the denominator of the integral in (10). Procedures 
are described in the literature^' '^ for computing the distribution of a 
finite quadratic form in Gaussian variables. Estimates of the error due to 
truncation can also be obtained rather readily. Unfortunately, the best 
such estimates obtained by the author showed that for small values of 
i9 or T, M must be taken quite large (50 or 60) to obtain answers guaran- 
teed accurate to two decimal places. Furthermore, the convergence of 
the computational schemes described^- ^ turned out to be very slow. The 
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Fig. 12 — Cut lines and contour in complex | plane. 
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following alteni.ative approach which can also be applied to finite sums 
of the form (5) was used to obtain the curves presented here. 

The function 11(1 — 'li^ajf' in the denominator of (10) has branch 
points at —ihj , where 

The ?/s are real positive (luantities, for the X's are eigenvalues of a real 
symmetrie positive definite kernel and must be real positive numbers. 
Line segment cut-lines are inserted in the complex J-phane from —ib-a to 
—ibij+i , J = 0, 1 , 2, ... as shown in Fig. 12. When y < 0, the value of 
(10) is zero as can be seen by closing the contour in the upper half plane. 
When y ^ 0, the contour of integration in (10) is displaced from the 
axis of reals to the contour, C, shown in Fig. 12. This dispUicement of 
contour is easily justified if 11(1 — 2i^ajf' is of exponential order less 
than luiity, a condition which wall be fulfilled in the examples to be 
treated. The change of viu'iable f ^ i^ rotates the contour of Fig. 12 by 
90° in the positive direction. If one now coUapses the closed contour 
curves about the cut-lines and takes proper care of the (convention al- 
ready set forth for the square root sign, there results, 



f(y) =E"^o(-i)V, 



where 



and where 



1 r'--"-'^ ^' c--"' dt 
^ 1,,, V\D{t) 1' 



A^- TTf^TTTl. ^ = 0,1,2,... 



D(t) =n(l -r)- (13) 



y=o 



bjl- 



D(t) is closely related to the Freeholm determinant (Reference 4, Chap- 
ter 11) of p. 

In the application to be treated below, 

D{t) = H{z) (14) 

where 

z = g{t) (15) 

is a non-negative monotone increasing real function of i for / ^ h) . De- 
note its inverse by i ^ h{z). Let 



178 THE BELL SYSTEM TECHNICAL JOURNAL, JANUARY 1958 



<^lc = 2(224+1 — 22t) 

for jfc = 0, 1, 2, . . . and let 

Z = Ck cos -KX -\- dk . 
Then straightforward substitution yields 



Ik 



Jo 



Y \ (z ~ Zik){z2k- 



Similarly, one obtains 



){Z2k+\ — z) 

f{y) =T.k^.{-\fh. . 



/c-0,],2,..., (16) 



(17) 



with 



/ 



Jo 



e-i>i^^'^h'(z) dx 



h(z) 



H{z) 



(Z — Z2k)(Z2k+l — Z) 



(18) 
(19) 



Equations (16) to (19) were used to compute the curves discussed in 
Section II. The denominators of the integrals in (16) and (19) have no 
zeros in the range of integration. By use of Gauss's method of numerical 
integration,^ evaluation of the integral at a; ^ and x = 1 where the 
denominator is an indeterminate form was avoided. In the applications 
made, it can be shown that for sufficiently large k, h and Jk decrease 
monotonely. Since the series (17) and (18) are alternating, an estimate 
of the error made by terminating the series at a finite value of k can 
be obtained. In all cases computed, it was never necessary to take k 
larger than 18, to obtain 1 per cent accuracy in the final result. 

V. DETERMINATION OF EIGENVALUES AND Hiz)* 

For stationary processes, the kernel of the integral equation (4) be- 
comes a difference kernel; that is, p(t, t') — p(t — t') where pix) is a 
positive definite function. The Fourier transform of p, namely 



w 



(/> - [ 

J— a 



2wifx 



{x) dx 



is non-negative and is the power density spectrum of the processes. 

* An alternative method of evaluating H{z) is described in Reference 12. 
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Analytic solutions to the integral equation (4) are known in this case 
only for a relatively small elass of kernels. Fortunately, this class is one 
of considerable interest in (communication applications. It is the class 
of p whose spectra w{f) are rational functions of P; i.e., ratios of poly- 
nominals in p. Such spectra are obtained by passing white noise through 
a finite passive physical electrical network with lumped constants. De- 
tails of the method of solution are given in References 9 and 10. It 
must be pointed out that, even in this case, solutions can be carried out 
practically only for polynomials of small degree. 

5.1 RC Noise 

If white Ciaussinn noise is apphed to a series RC circuit, the voltage 
across the capacitor has a power density spectrum proportional to 

where a = 1/RC is the nominal cut-off frequency of the circuit. The co- 
variance function corresponding to (20) is 

Pit) = fi-"'^'. (21) 

Solutions to (4) with this kernel are given in detail in both References 
9 and 10. 



Let 



Then 



/3 = ^ . (22) 



b, =-L[^^ + ./], ^-^ 0,1,2,..., (23) 

where the Zk are non-negative roots of either of the equations 

2 tan 2-/3 (24) 

zcoiz ^ -&. (25) 

If the zu are labelled so that 2o < 2i < 22 < ■ ■ - , then it is readily seen 
that Zk -- /.■(t/2), so that ht ^ />'^(tV8/3). The convergence exponent (see 
Reference U, ]). 14) of the secjuence h is therefore 5. It follows then (Ref- 
erence 11. 2.fi.5, p. 19) that D{t) as given by (13) is an entire function of 
order i. 

Now the func tion {c~'^/l3)[0 cos z + 2 sin z] Icos z + ^(sin z/z)], where 
z ^ -y/'Z^t — )3-, is an entire function of i of order ^. Its only zeros are 
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Eit the points t = h^ , k = 0, ], 2,. .. , given by (2:^), (24) iind (25). At 
i = 0, it hiLS the vtihie unity. It is, therefore, equal to D(l) as can be seen 
from Hadamard'a ractorization Theorem (Reference 11, 2.7.1, p. 22). 

The {]iiuntities necessary to evaluate (IG) and (19) are therefore all 
known for this case: 



H{z) = -— [/9 cos 2 — 2 sin z] 

Hz) = ^ t/5^ + ^'] 



COS 2 + ^ 



in zl 



and the zt are given by the positive roots of (24) and (25). 
The first two semi-in^'ariants of 7j arc found to be 

«] = Ey = 1, 

K2 = Ehj - If = / = 2. [4/3 - 1 + e-'\ 

5.2 RLC Noise 

If white Gaussian noise is applied to a series RLC circuit, the voltage 
across the resistor has a power density spectrum proportional to 



wq . . fUV , -y •,., (26) 



©'<-' 



u' + I — ] (w — Wo ) 



where w = 2jr/", Q = mL/R and m = 1/LC. Introducing parameters u 
and V defined by 



U' -\- V~ = Wo' 






uv = tou', Re u ^ 0, Re V ^ 0, 
one finds 

, , ^ 2(w + t')w^ 

and 

p(t) = -J- [ue-"^'' - vcr"^% (27) 

U ~ V 
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In the special case Q = h 

p{t) = (1 - coo 1 T IVo""" '^ (28) 

Solution of (4) with (27) or (28) as kernel is relatively straightforward 
by the methods of References 9 and 10, although quite laborious. Details 
can be found in Appendix 2. 

Suppose Q and wo are positive real quantities. Then the eigenvalues 
Xt = T/2bk are given by 

(,, ^i[r + ,,^ /. ^0,1,2, ... (2!)) 

2r 

where the z^ are non-negative roots of either 

. ,2 -A sin 2 ,. 2 I 2-. sin s/z' + s- ,..^^ 

2r cos z = {z - r) + (2 + r ) , (-iO) 

z Vz- + s- 



or 



Here 



sin -y/z^ + s^ 



2r cos 2 =- (2 — r) — — (2 + r ) • — . ■ (-iiJ 



The eigenfunctions belonging to roots of (30) are of the form 

Ak cos |(2a- + Vzk' + s'')t + B, cos K2t - Vzk- + s=)^ 
while those belonging to roots of (31) are of the form 

C\- sin M^A- + V^TTT^)^ + A- sin K^a - V'^7T~^')^ 

11 is interesting to note that when the X's :ire ordered in the usual way, 
the corresponding eigenfunctionsdonot in general alternate between even 
and odd functions of /. 

The infinite product (13) with the b's given by (29), (30) and (31) 
can be written in closed form by arguments similar to those used in 
Section 5.1. From (30) and (31), it is seen that asymptotically succes- 
sive zt are separated by ir/2, so that h grows like /-^ and one is again 
deahng witli an cntiie function of order L* For the pertinent ([uantities 

* More gt'iKMiilly, it fiiii be shown thiit for riitioiiiil spectrn if ii'ijj ^ /"-'' then 
X„ "^ H"^''. (Privtite commiiuioation to iiuthor liy A. Beurling.) 
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of (IG) and (19), one finds 



4r^ 



"(2) ^ -r^ U — r ) — 2r cos 2 

L 2 



+ (.^ + .yilL^^t?] (34) 

■ (2 ~ r ) 2r cos s — (s + r ) 

L 3 V2- + s^ 

/^(^) = ^^ [^' + r^] (35) 

with the Zk given as roots of (30) and (31). 
The first two semi-invariants of y are 

Wi = Ey = 1 

»^ = -' = ^ [2'- - 1 + «"" + 2 /^ sinh^ Vr"^^T^] . 
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APPENDIX 1 

Let N(i) = a cos W + b sin oiot, where a and h are independent Gaus- 
sian variates with mean zero and variance unity. Then y as defined by 
(1) is obtained by direct integration as 



y - aa + 0b^ 



where 



1 /, , sin s\ o 1 A sm sN 

2(1+^' '' = 2V-—)' '-''°^- 



Since y is the sum of independent random x variables, the density for y 
can be obtained as the convolution 



f(y) = f 

Jo 



V ~{xl2a) „-((y-i)/2|S) 

ax. 



'0 V'iirax V2Tr0(y - X 
The substitution x =^ (?//2)(l + cos d) in this integral leads to 
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Finally, if sin i^ ^ sin s/s, 

f(y) = secvw-^'^'VoCt"?/ tan ip sec v')- 

APPENDIX 2 

The power spectrum corresponding to the covariance (27) can be 
written as 

2(w + v)p^ 



w = — 



(p- — u^){p'^ — v^) 



where p = iu ^ 2inf. From Reference 9, then, solutions to (4) with 
the kernel (27) must satisfy the dilTerential equation 

or 



2 , .2 2,2 2(m + y) 



where 

2,3 2 2 4 

a ^ = uv = ilia • 

We choose a and (3 so that Re a ^ 0, Re /3 g 0. If a ?^ j3, then ^ is a 
linear combination of the elementary functions e"', e ° , e , e . 

It is easy to verify that if .^ is a solution to (4) with a kernel p{t, I') = 
p{\t ~ t' \), then <p{l) + ^(-0 and ^(0 - v(-i) are also solutions. 
We can, therefore, restrict attention to even and odd solutions of (4). 
On substituting 

^(0 = A cosh at -\- B cosh 0t 



into (4), one finds 

(ii) 



X aV 0'r 



T (u' - a=)(i'^ - «') (w- - /3=)(i'- - ^') ' 
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u cosh -^ + a sum -^ u eowli -^ -\- /Ssinh. -— 

A ^ , ^ + B ^-- ^ ^ 0, 

w^ — a^ w- — ,3^ 

(iii) 
V cosh -7^ + a sum -— v cosh -5- + /3 smh -— 
^ ^^ ^ + s ~^, — - 0. 

yj ~ a- V- — p~ 

The determinaut of the system (iii) must vanish. A bit of algebra shows 
this to be equivalent to 

2r cosh X + {x -r r ) + (x — r) — — = , (iv) 

where a; = (a + ^)(7'/2). It is not difficult to show that for positive wo 
and Q, this equation has roots only if a and /?, and hence x, are pure 
imaginary. Writing x = iz, (iv) become (30) and (li) yields (29). 
The substitution of 

ipit) = C sinh al -\- D sinh 0i 

into (4) again yields (ii) and equations analogous to (iii) with sinh and 
cosh interchanged. A similar analysis then gives (31). 

If a = /3, then from (i), ip must be of the form A cosh at + Bt 
sinh at or C sinh at + Dt cosh aL Substitution of these forms into (4) 
yields equations which cannot be satisfied for positive a)o and Q except 
by the trivial solution A =B=C = D = 0. 
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